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Abstract 

■ We consider the problem of the exact computation of the correlation functions of the 
j. I eight- vertex solid-on-solid model by means of the algebraic Bethe Ansatz. We compute the 

(SJ ■ scalar product between a Bethe eigenstate and an arbitrary state of Bethe type and show 

I that, in the cyclic case, it can be formulated as a single determinant of usual functions. It 

■ allows us to obtain determinant representations for finite-size form factors. By summing 
up over the form factors, we also give a multiple integral representation for a generating 

04 ■ function of the two-point function. 



1 Introduction 

In the domain of integrable systems, after the successful determination of the spectrum and the 
eigenstates of the Hamiltonian [101 123 EH U\ EH ES [2l] , a crucial problem is the computation 
of the form factors and correlation functions which are essential objects for the description of 
the full dynamical properties of the models. For truly interacting integrable models (i.e. for 
models which are not equivalent to free fermions) , this problem represents a real challenge and 
has been intensively studied during the last decades (see e.g. [ST, 19, 3?. IM l [HI ] 1^ W5[ Hil [35 1 
[TiEIllZaiiaiMllMllMlSQlEaiiaiiailllllZllia^ several methods 

have been developed so far, and important progresses have been made recently, in particular 
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for the archetypal integrable models such as the isotropic (XXX) and partially anisotropic 
(XXZ) Heisenberg chains, or the one-dimensional quantum Bose gas. 

Among these methods, one [62l [63], based on algebraic Bethe Ansatz (ABA) [251 l23] . 
has led to important recent developments [MIEZllMIMlEaEaiMlEaEZllMlEHlEZ]. It 
relies on the solution of the quantum inverse problem [621 [771 [39] and on a compact and 
explicit expression, in the form of a determinant of usual functions, for the scalar products 
between a Bethe eigenstate and an arbitrary state of Bethe type [90] . The combination of 
these two results leads to simple determinant representations for the form factors (i.e. for 
the matrix elements of local operators in the basis of eigenstates of the transfer matrix) in 
finite volume [62j. Although the exact representations which are obtained for the physical 
correlation functions (and in particular for the two-point functions) have a less simple form 
[Ml [HU [601 [59l [S] 5 it is nevertheless possible to obtain quantitative results for these quantities 
as well, for instance by summation over the corresponding form factors series. This can 
be achieved either by numerical methods relying on the efficient, explicit aforementioned 
determinant representation for the form factors |11[ [2U1 or, in the asymptotic regime, from 
the analytic study of the series [551 EH EHl [Ml [57] • 

However, until now, this method has been essentially developed in the relatively simple 
case of the periodic XXZ spin-1/2 Heisenberg chain or some of its variants (open chain [521 153] . 
higher spins [SH [JOl [22] . . . ) , or in the even simpler case of the Bose gas model [Ml [671 EZ] ■ 
In fact, for more complicated models, the solution of one of the two aforementioned basic 
ingredients (resolution of the inverse problem and determinant representation for the scalar 
products of Bethe states) is often missing. 

In particular, a real challenge would be to adapt this method so as to find explicit and 
manageable expressions for the correlation functions of the completely anisotropic (XYZ) 
Heisenberg chain. The latter is a natural generalization of the XXZ Heisenberg chain, and is 
related to the eight-vertex model of statistical physics (which in its turn is a natural general- 
ization of the six- vertex model related to the XXZ chain). It was solved for the first time by 
Baxter [3l [H [5l [6] (see also [59] for the computation of the correlation length, and [92l [28] for 
the ABA formulation). In fact, it is worth mentioning here that, due to the non-conservation 
of the spin, the XYZ (or eight-vertex) model is not directly solvable by Bethe Ansatz. Nev- 
ertheless, Baxter managed to construct the eigenstates of this model by relating them (via 
the so-called Vertex-IRF transformation) to the Bethe eigenstates of another model of solid- 
on-solid (SOS) type. The latter, which is often called eight- vertex SOS model (8VS0S), is 
a L-state IRF (interaction-round-a-face) model on a square lattice and is solvable by Bethe 
Ansatz [U1[2H]. However, due to the need of the Vertex-IRF transformation to come back to 
the XYZ (or eight-vertex) model, the resolution process is therefore much more complicated 
than in the XXZ case, and the exact computation of the correlation functions is in this case 
still a widely open problem. In particular, through the aforementioned method based on 
ABA, no significant results in this direction have been obtained so far (see nevertheless, by 
means of another approach using (7- vertex operators [53], some first results based on those of 
[72] in [TOl [69], as well as an attempt of a direct approach in [57]). 

Actually, the combinatorial complexity issued from the use of the Vertex-IRF transfor- 
mation is not the only problem for the computation of the XYZ correlation functions. In 
fact, even in the simpler case of the related SOS models [681 ESl [82], the ABA approach 
to correlation functions presents some difficulties, in particular concerning the obtention of 
a compact and manageable formula for the scalar products of states. Let us recall at this 
point that, in the XXZ case, the scalar product between a Bethe eigenstate and an arbitrary 
state of Bethe type can be written as a determinant of usual functions [90J, and that the ABA 
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approach to correlation functions relies strongly on such a representation [62]. The latter orig- 
inates (see for instance the derivation of |62) ) from the explicit determinant representation of 
the partition function of the six- vertex model with special boundary conditions, the so-called 
domain wall boundary conditions [641 142j . The main problem in the 8VS0S case is that the 
analogous partition function with domain wall boundary conditions does not seem to admit 
a representation in the form of a single determinant [84:\ [HI] (see also [33] ) . This is due to the 
fact that the i?-matrix of the model depends (compared to the simpler six-vertex i?-matrix) 
on an extra dynamical parameter related to the fluctuation variable (height) of the model. 
This i2-matrix hence satisfies the dynamical version of the Yang-Baxter equation [26\ \29\ [28] 
instead of the usual one. As a result, the dynamical parameter is subject to some shifts, which 
prevents one from reducing the sum of determinants obtained in [M] for the 8VS0S partition 
function to a unique one as in (see nevertheless [31[ [30] for a model with a reflecting 
end, for which the partition function indeed reduces to a unique determinant of Izergin type). 
Hence, a mere generalization of the process described in |62j seems to be not so easy. 

This is the problem that we tackle in this article, namely to set the bases of the ABA 
approach to correlation functions in the dynamical case, i.e. in the case of the periodic 8VS0S 
model. Building on Rosengren's representation [84| of the partition function with domain wall 
boundary conditions as a sum of determinants, we obtain a similar representation (i.e. as a 
sum of determinants) for the partial scalar product, evaluated at a fixed value of the height s, 
of a Bethe eigenstate with an arbitrary state. Hence, the fact that we do not a priori obtain a 
single determinant may be a problem for the computation of correlation functions. However 
we show that, in the cyclic case (i.e. when the parameter rj of the model is rational and the 
space of states is finite-dimensional), the true scalar product of Bethe states, which is obtained 
by summing up over all values of the dynamical parameter, can in fact be represented as a 
unique determinant. In the same way, the finite-size form factors (i.e. the matrix elements, 
in the basis of eigenstates of the finite-size transfer matrix, of local operators labeling the 
difference of height between two neighboring vertices) can be written as a unique determinant 
which has a similar form as in XXZ. This opens the way to the computation of correlation 
functions, for instance by summation over the form factor series. In particular, we define a 
generating function of the two-point function and show that, in the cyclic dynamical model in 
finite volume, the latter admits a multiple integral representation very similar to the master 
equation representation that was obtained in [60j and used in [55] to derive the long-distance 
asymptotic behavior for the two-point function in the XXZ case. 

The content of the article is the following. In Section [2] we define the 8VS0S model and 
recall the main steps of its algebraic Bethe Ansatz resolution. In Section [21 we compute the 
scalar product between a Bethe eigenstate and an arbitrary state of Bethe type. In Section [H 
we solve the quantum inverse problem in the dynamical case, namely we express local spin 
operators in terms of elements of the dynamical monodromy matrix. In Section [Sj we compute 
the finite-size form factors. In Section [6l we formally sum up the form factor series and obtain 
a multiple integral representation (master equation representation) for a generating function 
of the two-point function of the (finite-size) cyclic dynamical model. Details and technicalities 
are gathered in a set of appendices. 

2 The 8VSOS model and algebraic Bethe Ansatz 

Let us consider a two-dimensional square lattice with N x M elementary square cells (or 
faces), with periodic (toroidal) boundary conditions. A height s is attached to each vertex 
(site) of the lattice, so that heights on adjacent sites are restricted to differ by ±1, and a 
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complex parameter Ui (resp. ^j) is attached to each column i (resp. line j) of cells. To each 
height configuration (si, S2, S3, S4) (with \si — S2I = |s2 — s^l = |s3 — S4I = |s4 — si| = 1) 
around an elementary face of the lattice is associated a statistical weight W^'l,^^ which can 
be understood as a matrix element of some matrix R depending on the difference of the two 
corresponding parameters: 

s + e' 







J. 

Ui 



s + ei + Cj 
= s + e[ + e'^ 



Here , , , e^- € {+1,-1}, and the matrix R{u; s) G End(y (g> V), where F ~ C is a 
two-dimensional vector space with basis (e+ , e_ ) , is such that 

i?(n;s)(e,/ ®e,/) = Riu; s)l}'2 {^e^ ^ e,,) . (2.1) 



E 

<:i+<:2=<:'i- 



The i?-matrix of the SOS model admits only six non-zero elements, which can be param- 
eterized as follows: 



R{u] s) 



with 



b{u; s) 
b{u; s) 



/I 0\ 

b{u; s) c{u; s) 
c(n; s) b{u; s) 

Vo 01/ 

[^ + 1] M 
[s] [u + iy 

[s - 1] M 

[s] [u + 1] 



(2.2) 



b{u; -s), 



c[u; s 



c[u; s 



[s + n] [1] 

[s] [u + iy 

- [1] 
[s] [u + l] 



c(u; 



(2.3) 
(2.4) 



The parameters u and s are respectively called spectral and dynamical parameters. The 
function u 1— )• [u] is an entire, odd and quasi-periodic function of quasi-periods 1/r] and 
t/t], where rj and r are two parameters such that Q'r > (see Appendix [X] for the precise 
definition of this function). Graphically, the above non-zero matrix elements correspond to 
the six following statistical weights: 



s + s + l 
+ + 



s+l + s+2 
a{u; s) = 1 



s + s + l 



s-1 + s 
b{u; s) 



s - 1 

+ 



s-1 + s 
c{u; s) 



s-1 



s-1 - s-2 
a{u; s) = 1 



s 
+ 
s + l 



s - 1 

+ 

s 



b{u; s) 



s 
+ 

s + l - s 



s + l 



c{u; s) 
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The i?-matrix ()2.2p satisfies the dynamical Yang-Baxter equation [36\ [26] (which is equiv- 
alent to Baxter's star-triangle relation for the Boltzmann weights W) : 

Ruiui -U2;s + h-i) Ruiui - U3; s) R23{u2 - U3; s + hi) 

= R2z{u2 - U3; s) Rl3{ui - Us; S + /l2) Ri2{ui - U2; s), (2.5) 

with 

-0- <^-^> 

This equation should be understood as an identity for meromorphic functions of ui,U2,U3, s 
with values in End(y i^V^V). The indices label as usual the space of the tensor product 
on which the corresponding operator acts. The i?-matrix (|2.2p also satisfies the following 
properties: 

• zero weight: [Ri2{u; s) , hi + /12] = 0, 

• unitarity: Ri2{u; s) R2i{—u; s) = Id, 

• crossing symmetry: a\Ri2{—u — 1; s — hi)a\ — ^^^J" ~ ^21 (^i 

in which, in the expression Ri2{—u — 1; s — /ii), the hi operator should be understood 
as acting to the right of any other operator involved in the definition of the i?-matrix. 

The transfer matrix of the model corresponds to the product of all statistical weights along 
a column of elementary cells of the lattice. A given allowed configuration of heights along a 
vertical line of vertices of the lattice (a state on which the transfer matrix acts) corresponds to 
a A^-tuple of heights (si, S2, • • • , sat) such that Sj+i — Sj = ±1, i = 1, . . . N (with the convention 
sn+i = Si), i-e. to a configuration (5,6^^, . . . ,ee^) with s = si and = Sj+i — Si such that 
ei + • • • + eAr = 0. Hence, the space of states of the model can be seen as the space of functions 
Fun('H[0]) of one complex variable (the height s) with values in the zero- weight space ?^[0], 
where H = V®^ . 

In the general (unrestricted) SOS model, r/ is arbitrary and the dynamical parameter s 
may take an infinite discrete set of values, i.e. belongs to some set Csq = sq + Z for some 
arbitrary parameter sq- Hence, in that case, the space of states is infinite-dimensional even for 
a finite lattice. Usually, however, one considers situations for which the height s is only allowed 
to take a finite set of values and the space of states is finite-dimensional. This happens for 
rational values of ru i.e. for rj = r/L with r, L being relatively prime integers such that L > 
(or more generalljo, as in [6], when there exist two integers ri and r2 such that Lr] = ri +r2r). 
In that case, that will be referred to as the cyclic case (or cyclic SOS (CSOS) model [83j), the 
statistical weights of the model (i.e. the elements (|2.3p - ()2.4p of the i?-matrix) are periodic of 
period L, and the space of states H^^^ corresponds to the space of functions ijj : Csg 1-^ 'H[0] 
such that ip{s + L) = il){s) or, equivalently, to the space of functions : C^^, 1— )• ^[0], with 
=so + Z/LZ. 

The algebraic Bethe Ansatz approach to the SOS model, which enables one to diagonalize 
the transfer matrix of the model, has been developed in ^28j, based on the study of repre- 
sentations of Felder's dynamical quantum group E-r^rjish) associated to the dynamical 

'^This situation is actually equivalent to the case rj = r/L with r being the greatest common divisor of ri 
and r2 (see Appendix [X]). 
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i?-matrix ()2.2p . The central object of this approach is the monodromy matrix, defined as the 
following ordered product of i?-matrices along a column of elementary cells of the lattice: 



Ta,l...Niu;(,l, ... ,(,n;s) = RaN{u - ^N] S + hi -\ h hjy-l) ■ ■ ■ Ral{u - Cl] S 

^A{u;s) B{u;s) 



It acts on Va T-L, where ~ is another copy of V usually called auxiliary space. In 
this framework, the entries A, B, C, D of the monodromy matrix are linear operators acting 
on T-L, and the commutation relations of these operators are given by the following quadratic 
equation on (8) ® 

Raiaziui - U2; S + hi,„N) Ta^,l...N{ui; s) Ta2,l...N {U2] S + ha^) 

= Ta2,l...N{u2;s) Ta^,l...N{ui;S + ha2) Raiaii^l -U2;s), (2.8) 

which is a consequence of the Yang-Baxter relation (j2.5p . Here hi,,,M = hi + . . . + hjsf. 

It may be convenient to define, from the above monodromy matrix T G End(VIi (8) Ti), the 
following matrix of finite-difference operators: 

nn)=fe bIII =^(^;^)fn GEnd(T4®Fun(?^)), (2.9) 

where r<j s = {s + 1) Tg, and the action of s and Tg on functions / G Fun(?^) are given as 

[sf]{s) = sf{s), [%f]{s) = f{s + l). (2.10) 

This defines an operator algebra whose commutation relations follow from (|2.8|) (see |29' '28]). 
In this picture, the (operator) transfer matrix is t{u) = A{u) + D{u). It is easy to see [29^ i28j 
that these transfer matrices preserve the space Fun(?^[0]) of functions with values in the zero 
weight space ?^[0] (the subspace of H on which hi,,,N vanishes), and that they commute 
pairwise on Fun(?^[0]): ^{u),t{v)] = on Fun(-?^[0]). 

The structure of the i2-matrix implies that there exists a reference state | ) = e+ (8) ■ ■ ■ (8) 
e_|_ G 7^ which is an eigenstate of A{u] s) and D{u; s). More precisely, 

A{u- s)| ) = a(n)| ), D{u; s)\Q) = d(n)| 0), (2.11) 

with, in our normalization. 

Starting from this reference state, one can construct states in the zero- weight space Fun('H[0]) 
as 

s ^ ip{s)B{vi; s)B{v2; s - I) . . . B{vn; s - n + 1)| ), (2.13) 

where vi, . . . ,Vn are arbitrary spectral parameters, (p is an arbitrary numerical function of the 
dynamical parameter and n is such that N = 2n + i^L for some integer K (or = 2n in the 
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unrestricted case). Bethe states correspond to states (|2.13p for which the function ip takes 
the form 



v'-(«)='^^n A' (2.14) 

depending on some complex parameter oj. Such states will be denoted as | {^i, . . . , u„}, w ), 
or simply as | {v}, u ) when there is no ambiguity concerning the set {v} = {vi , ■ ■ ■ , f «}• Note 
that, whereas u can a priori take any complex value in the unrestricted case, it is no longer 
the case in the cyclic case: for such a state | {v},uj) to belong to H^^^ (i.e. to be L-periodic 
as a function of s), one should impose a; to be such that (— l)'""a;^ = 1. 

It can easily be shown, using the commutation relations coming from ( |2.8p . that, if 
vi, . . . ,Vn are off-diagonajl solutions to the following system of Bethe equationql 

then 

t{u) I {v},uj) = t(u;{v},uj) I {v},uj), (2.16) 

with 

r(n; {.}, u) = a(n) fl + JJ t^i+A, (2.17) 

^^■^ [vi - u\ [u - Vl\ 

Similarly, setting (0| = |0)^ and 

^.(s)=a;-n^' (2.18) 
i=o ^ J 

one can define Bethe states ( {v},uj \ in the dual space of states as 

{{v},uj\ : {0\C{vn;s-n)...C{v2;s-2)C{vv,s-l)^Us)- (2.19) 

Then, if the set of spectral parameters {v} = {vi, . . . is an off-diagonal solution of the 
system of Bethe equations ()2.15p , one has 

{{v},oj\tiu) = Tiu;{v},u){{v},uj\, (2.20) 

with r(u; {v},ijj) given by (j2.17p . 



3 Scalar product of Bethe states 

The next step towards the calculation of form factors and correlation functions, after the 
determination of the eigenvectors of the transfer matrix, is to obtain manageable and compact 
formulas for their scalar products. To compute these scalar products we use here the method 
proposed in [62], which is based on some recursion relation initiated from the determinant 

''i.e. such that, for all i < j, i]Vi 7^ r/Vj mod Z + rZ. 

^These equations are also valid in the unrestricted case for N = 2n with the convention rK = 0. 
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representation of the partition function with domain wah boundary conditions |42j . This 
recursion uses the expression of the operators B and C in the so-cahed F-basis, a basis of the 
space of states in which these operators become quasi-local [76j. 

In the case of the dynamical SOS model, the F-basis has been exphcitly obtained in [T], 
whereas an expression (unfortunately not as a single determinant but as a sum of determi- 
nants) for the partition function of the model with domain wall boundary conditions has been 
obtained in [811 184j (see Appendix [Bj. The fact that one does not know a compact formula 
(as a single determinant) for the partition function is a difficulty for the computation of the 
scalar product, but the approach of [62] can nevertheless be achieved and, at least in the 
cyclic case for which the space of states is finite-dimensional, the result can be presented as a 
single determinant. 

Let ( {u}, I be a dual Bethe eigenstate given in terms of a function (^(^^ of the fom (|2.18p 
with complex parameter w^, {u} = {ui, . . . ,n„} being solution of the corresponding Bethe 
equations, and | {v},uj^ ) G H^^ be a state of the form ()2.13p given in terms of a function ip^^^ 
(|2.14p with complex parameter uj^, {v} = {vi, . . . , Vn} being a set of arbitrary parameters. In 
the cyclic case, we consider the following scalar product between these two states: 

S„({n}, uju] {v}, ujy) = {{u},uju\ {v}, ) 

1 „ _ 

Vuuis) VujAs) Sn{{u}; {v}; s) 



L ^ 

"so 



= 1 E ^u'<i]}^^^^Sn{{u\Av\.s\ (3.1) 

where S'„({'u}; {u}; s) is the partial scalar product at a given height s defined as 

5„(M; {v}; s) = ( |C(u„; s - n) . . . C(mi; s - \)B{vv. s) . . . B{vn; s-n + l)\0). (3.2) 

Following the method proposed in [62], ()3.2p can be calculated in the F-basis because the 
reference state and the dual reference state are left invariant under the action of the F-matrix 
which induces the corresponding change of basis (see \76\ l62l |1] for more details about the 
F-basis). (|3.2p is therefore equal to 

Sn{{u}; {v}; s) = (0 |C'(u„; s - n) . . . C(ui; s - l)B{vi; s) . . . B{vn; s-n + l)\0), (3.3) 

where C and B stand for the expressions of the corresponding operators C and B in the 
F-basis (see Jjlj for details): 

Bin; s) = f " ^r'U i^ . m 




+ hl...N] ' [s + l:l^^hl][u-ii + l]3^ 

ir[ \s\[u-ii + i\ j^i y i/j^.j 

As already mentioned, (|3.3p can be computed by recursion. For k G {0, 1, . . . n}, we 
consider the following quantity: 

g\%,,.../„{{u}; {vi, . . . ,vk}; s) 

= (0 |C(u„; s-n)... C{ui; s - l)B{vi-s) . . . B{vk; s-k + 1)| 4+i, • ■ • ,4, ), (3.6) 
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where | ik+i, ■ ■ ■ ,in) denotes the state with N — k down spins at sites ik+i, ■ ■ ■ ,^n- For k = n, 
the quantity G^^\{u}; {vi, . . . ,Vn}; s) corresponds to the partial scalar product (|3.3|) . For 
k = 0, G^^^ ^^({ii};0;s) can easily be computed, by means of the expression ()3.5p of the 
operator C in the i^-basis, in terms of the partition function with domain wall boundary 
conditions (|B.ip on a lattice of size n x n (we also use the fact that the F-matrix leaves 
invariant the state | ) so that (jB.ip can be written directly in the F-basis) : 

0; = ( |Ci...7vK; s - n) . . . Ci...N{ur,s - l)\h, . . . Jn) 



X Z„(ni,...,u„;^£,,...,^^^;s-n). (3.7) 
The functions (|3.6p admit the following recursion relation: 

GfJ+^,...,inii^y^ • • • , vk}; s)= Gf~\{{u}; {vi,..., Vk-i}; s) 

45^4+1, -A 

X {(k, ... ,in\B{vk;s - k + l)\ek+i, ■■■ Jn), (3.8) 
with, using the expression (|3.4p of the operator B in the F-basis, 

[l][s + N -2n + k + Vk- ^ij 



{ik,---,L\B{vk;s-k + l)\ik+ 



[s + N-2n + k] [vk - e^J 



ii _ 1 _ ^„ 1 • (^-9^ 



The recursion relation (j3.8p with initial condition (|3.7p is explicitly solved in Appendix O 
Here we merely give the result for the partial scalar product ()3.2p . 



Proposition 3.1. For N = 2n + KL with ^^ integer, let {ui, . . . , Un} be a solution of the Bethe 
equations (j2.15p uiii/i complex parameter oju, and {f i, . . . , Vn} be a set of arbitrary parameters. 
Then the partial scalar product S„({n}; {u}; s) can be represented as the following sum of 
determinants: 



s..(M;M;-=) = ,i„,, i'_7„i' ,,,, n 



s-m U\s-j\ n?=i 



X 



[7]"[|n| - 1^1 + 7 + s] [s + j] UjKki'^j - Uk][vk - Vj] 
5,5c{l,...,n} j^S ^ ^ ^' t=l ) j(zs I t=l 



^ [j + s-\S\ +_\S\] . _ . 

[s-\S\ + \S\] « ^ ^ ^ 

/n i/iis expression, 7 is an arbitrary complex parameter (the result does not depend on 'j), 
\u\ = ui + ■ ■ ■ + Un, \v\ = fi + • • • + Vn, and the elements of the n x n matrix are given by 
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The sum in ()3.10p runs over all subsets S and S of {1, . . . , n}, l^l and \S\ being the cardinality 
of these subsets, and {v — 6^^} = {vi — 6f ^ , . . . ,Vn — ^n^}, with 



6: 



ss 



1 if k ^ S and k ^ S, 

-1 if k^ S and ke S, (3.12) 
ifk^SuSorkGSnS. 



Remark 3.1. The formula ()3.10p is also valid for any generic SOS model (i.e. for generic r/) 
in the case N = 2n (U = 0). 

Remark 3.2. The apparent poles at Vk = Uj are removable due to the fact that {u} satisfies 
the Bethe equations. The apparent pole at |n| — + 7 + s = is also removable due to 
Remark IB. II 

In the case when r] is rational, i.e. if there exists a positive integer L such that Lr] is 
integer, the sum of determinants in (|3.10p can be reduced, by a similar argument as for the 
partition function in [84) . to a sum over only L terms. Indeed, using ()A.9p and the periodicity 

of the theta-function to re-express the ratio ^'^,''''^1 Ji^^^'jf^^ in (13.101) . one obtains the following 
result: 

Corollary 3.1. Suppose that there exists a positive integer L such that Lr] is integer. Then, 
with the same hypothesis and notations as in Proposition \3.1l the partial scalar product 
Sn{{u}; {v}; s) can be written as a sum of only L terms: 

^AW, in, s) ^1^1 -\v\+j + s] ii [s + j- 1] n,<,K - nk][vk - V,] 

X ff^-^^^^^ ^et K)(W,.„; W)], (3.13) 

with 

n 




-^d{v,)ll[ut-vj-l]. (3.14) 



t=i 

Here we have set q = e^'^*'' and [u]^ = Oiirju^Lr). 

Although we do not obtain a single determinant, the representation (j3.13p seems never- 
theless more convenient than ()3.10p for the computation of the correlation functions, its main 
advantage being that the number of terms remains finite in the thermodynamic limit. 

At this point, we would like to make the following remark. As we have seen, we did 
not succeed to represent the partial scalar product ()3.2p as a single determinant as in |90j . 
However, it is worth stressing that the state 

B{vi]8)B{v2;s-l)...B{vn;s-n + l)\Q) G ?^[0] (3.15) 

is not a Bethe state: even when the parameters {f } satisfy the Bethe equations, the shifts 
of the dynamical parameter that appear when commuting A{u\ s) or D[u; s) with B{u; s it 1) 
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prevent this state from being an eigenstate of, for instance, A{u; s) + D(u; s). Therefore, it is 
not very surprising that the occurrence of dynamical shifts also prevent us from re-expressing 
the quantity ()3.10p or ()3.13p (and the partition function (jB.ip as a particular case) as a single 
determinant as in the non-dynamical case. In fact, if we consider instead the true Bethe state, 
i.e. the function (|2.13p . we shall see that the natural scalar product of functions (|3.ip can, in 
the cyclic case, and for one of the set of parameters being solution of the Bethe equations, be 
represented as a single determinant. 

Let us therefore now consider the scalar product ()3.ip . Using the L-periodicity and the 
fact that the right hand side of (j3.10p does not depend on 7 and hence that a convenient choice 
can be made for this parameter, we can extract the dependance in the dynamical parameter s 
from the sum of determinants and factorize the latter as a single determinant. More precisely, 
we get from ()3.10p 



= bF n,<.K -dN 



X 



n ni". - + u n ni-. - ^= - n 



L f-^L [\u\ - \v\ + 7 + s] [s-\S\ + \S\] 



X det [J\f^{{u};{v -5^^})]. (3.16) 



Setting 7 = — |n| + and performing a change of indices in the sum over s, we obtain 

[7 + s]l 1 Ut=id(.ut)d{vt) 



Sn{{u},UJu;{v},UJy) = 1^ ^ 



X n n[^* - - 1] \ W,{{uh {v - 5^''})] • (3.17) 



t=l 



Using the linearity of the determinant, one can now express the last sum in (|3.17p as a single 
determinant. We have the following result: 

Theorem 3.1. For N = 2n + with integer, let {ui, . . . ,n„} be a solution of the Bethe 
equations (|2.15p with complex parameter cou, and {vi, . . . , Vn} be a set of arbitrary parameters. 
Then, the scalar product S„({u}, w^; {u}, w^) (|3.ip between the Bethe eigenstate {{u},uju \ o-nd 
the state \ {v},uji, ) G H^^^ of the form (j2.13p . (j2.14p with parameter can be represented as: 

o .r 1 r 1 \ J 1 '^«[7 + 'S]l Ilt=id{ut) 

I sfcL ^« J nj<fcK- - uk][vk - vj] 



X det [n^{{u},u;u; {v},^^v)] , (3.18) 
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with 7 = —\u\ + \v\ and ^l-y given by 



{-ly^ J [Uj - Vj + 7] [Uj -Vj 

[7] \ [Ui - Vj] 

1 J [Uj - Vj + 7] UJu [Uj -Vj+'J-l] 



/r 1 r 1 nT " -L lUi — Vj + 7 tij — Vj + ^ + L\\ , , -T-r , 



[7]\ [n,-t;,] a;„ \n,-v,-\\ j" ^ '^fj^ ^ J ^ ^ 

From this determinant formula, it can be shown (see Appendix[D|) that two different Bethe 
eigenstates are orthogonal. On the contrary, taking the limit in which {v\ = {u} with = lVu 
in (|3.18p . we obtain the formula for the "square of the norm" of a Bethe eigenstate (in that 
case 7 = 0): 

o .r 1 r 1 \ Ut=id{ut) ■detn[no{{u},uJu;{u},uJu)] 



Uj^k[uj-Uk] 



det [HM)]: 



Uj^k[uj-Uk] 

(3.20) 



with 



[$(M)] = dj^^og' ^{uj) + J2Kiu, -ut)^- k{u, - Uk). (3.21) 

Here we have defined the even function 

^(„) = fe^_^. (3.22) 

[u - IJ [u + IJ 



Note that the formula ()3.20p is very similar to its XXZ analog \35\ [M] . 

Remark 3.3. In the case ^^ = (n = N/2), the above formula does not explicitly depend on 
the periodicity L of the model. Hence, one can easily take the limit L — ?> oo, extending by 
continuity the validity of ()3.20p to all values of rj, including irrational ones. In that case, the 
cyclic condition has only been used so as to avoid the subtleties of dealing with an infinite- 
dimensional space of states. 



4 Solution of the quantum inverse problem 

In the ABA framework, the Bethe states are constructed by a repeated action of n non-local 
operators B on the reference state |0). In order to compute form factors and correlation 
functions, one needs to be able to act explicitly with local operators on such states, which 
may seem not so easy (at least by direct computation) due to the mixture of local and non- 
local operators in a same expression. A way to solve this problem was proposed in [62]: it 
relies on the solution of the quantum inverse problem, which consists in expressing the local 
operator we consider in terms of the generators of the Yang-Baxter algebra; as a result, it 
is possible to compute the action of this local operator on a Bethe state by using only the 
quadratic Yang-Baxter commutation relations given by the ii-matrix. 
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The method described in |62^ I77j . which rehes on the fundamental property that the 
i?-matrix evaluated at coincides with a permutation operator, can also be applied to the 
dynamical case. Indeed, the structure of the i?-matrix ()2.2p implies that, for any value of s, 



R^jiO;s)=Pij, (4.1) 

where Pij is the permutation operator of spaces i and j. The explicit reconstruction of local 
operators in the dynamical SOS case involves however some subtleties (with respect to the 
simplest XXZ case) due to the presence of the dynamical parameter s. 

Let E'j'^ be the elementary matrix, acting on the i-th space of the tensor product T-L = 
, with elements {E'^^)jk = SfS^, where a and /3 are equal to ±1. In order to adapt the 
method of EZ] to the dynamical case, we first note that 

h,Ef = Ef{h, + a-l3). (4.2) 

We shall also use the following lemmas: 

Lemma 4.1. For any value of s, we have the identity 

Tai,l...N{Ci;s) TaA...N{u] S + ka^) = Ta,2...NaA'^; S + hi) T^^ ,l...Ar (^1 ; S + /la), (4.3) 

in which Ta^i...N (respectively Ta^2...Nai) is the monodromy matrix of a chain of N sites labelled 
(in this order) by 1,2,...N (respectively by 2,3, . . . N,ai) with inhomogeneity parameters 
Ci, 6, • • • (respectively ^2, Cs, • • • ^N,Ci)- 

Proof — Applying the quadratic commutation relation (j2.8|) on the l.h.s. of (|4.3|) . we get 

Tai,l...N{(i;s)Ta,l,„N{u;S + hai) = Ra^ai^l - U; S + hi,,,^) 

X Ta,l...N{u; s) Ta^,l...Ar(6; S + ha) RaiaiCl -U;s). (4.4) 

Then, using the fact that i?aii(0) = Paii and that R is of weight 0, we have 

Taul...N{^i;S + ha) Raiai^l - U] s) = i?la(6 " U] s) T^j ,l...Ar (6 ; S + ha). (4.5) 

On the other hand, the i?-matrix being unitary, we have 

-Raia(6 -U;S + hi,„N)Ta,l...N{u; s) Rla{Cl - U] s) = Ta,2...NaAu; S + hi), (4.6) 

which ends the proof. □ 

Lemma 4.2. For any value of s, we have the following identity between products of mon- 
odromy matrices: 

Tai,l...N{Ci;s) Ta2,l...N{£.2; S + ha^) . . . Ta^,l...N{^u S + /Iq^ + /Iqj H h /iai_i) 

= Tai,i...Naia2...a,^i {£,u S + /li + /l2 + • • • + /ij-l) 

X 7ai_i,i-l...Araia2...ai_2(?«-i; S + hi + h2 + ■ ■ ■ + /li-2 + haj ■ ■ ■ 

■■■Ta2,2...NaAi'i\S + hi + /lag H V ha^)Ta^,l...N{il] S + ha^ ^ ^^aj, (4.7) 



where the notations are similar to Lemma \4.1\ (i.e. for instance Ta^^i...Naia2...ai-i denotes the 
monodromy matrix of a chain of N sites labelled by i, . . . N,ai,a2, ■ ■ ■ fflj-i with inhomogeneity 
parameters ^i, ^j+i, . . . Cat, 6, 6, • • • ^i-i)- 

Proof — (|4.7p can be proven by induction on i using (|4.3p . □ 
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The solution of the quantum inverse problem, i.e. the reconstruction of local operators 
E"^^ in terms of the entries of the monodromy matrix (|2.9|) . can be formulated as follows: 

Theorem 4.1. A local operator E"^^ acting on the i-th space of the tensor product % = 

can he expressed in terms of the entries of the monodromy matrix ()2.9p in the following way: 

i-l 1 

Ef = n ^fe) • ^paiii) ■ n Wk)] • (4-8) 

k=l k=i 

Hence, a product of local operators on adjacent sites admits the following reconstruction: 

fc=l k=i k=i+j 

Proof — Let us first prove ()4.8p for i = 1. Writing T^o(^i) as a trace over an auxiliary 
space, and expressing the monodromy matrix (|2.9|) in terms of products of i?-matrices and 
shift operator r^, we get 

fpaiil) = tra[f,,i...jv(ei;ei, • • ■,iN)Ef] 

= iTa[RaN{il - ^N^S + hi + ■ ■ ■ + h^-l) . . . -Ra2(6 " 6; S + h) Pal Ef], 

in which we have used (j4.ip . Passing in this expression the operator Ef from right to left 
by means of ()4.2p . we obtain 



TfSai^i) = tra[Ef RaN{£.i - ^n;s + hi -\ h hj^-l + a - /3) . . . 

...Ra2{Cl-C2;S + hi + a-/3) Pal f>+"-^] 
= Ef tTa[Ta,l...N{Ci;Cl, . . . « " /?) 

= Ef?^-n{(i), (4.10) 

which ends the proof of ()4.8p in the case i = 1. 

In the general case, one proceeds similarly by writing 131=1 ^i^k) ' Tf^ai^i) as a trace over 
auxiliary spaces: 

i-l 



k=l 

= trai...ai_iaJrai(Ci;s) • . . Ta^,^{Ci-i;S + ha^ -\ h /la,_2^ 



X TaM;S + ha,+--- + /ia._Jf.'-+-+'--+'- Ef]. (4.11) 

Passing the elementary matrix Ef from the right to the left, first through the shift operator 
% ^ g^j-^j then through the product of monodromy matrices, one obtains 

i-l 

■ Tp^m = Ef tVa,...a,_,aATaAil-^S + a - fi) . . . 

k=l 

■ ■ ■ Ta,_i{Ci-l]S + ha^-l h /la,_2 + a - f3) 

X TaM^+ V + • • • + ha^_, + a - /3) ^^.,+-+h.,_,+h.,+a-P^^ ^^_^2) 
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Here we have used ()4.2p . and the fact that the product of monodromy matrices could be 
expressed in a convenient way by means of Lemma 14.21 so that its commutation with Ea^ 
(which becomes i?"^ by this process) can be performed similarly as in the case i = 1. Finally, 
it remains to note that ()4.12p can be rewritten as 



i-1 



n ■ fp^m = Ef n ^(^'«)' (4.13) 



k=l k=l 



which ends the proof of ()4.8p . 

()4.9p can be proven by induction from ()4.13p . □ 



5 Determinant representation for finite-size form factors 

In this section we consider finite-size form factors of the 8VS0S model in the cyclic case, i.e. 
matrix elements of local operators between Bethe eigenstates of the finite-size transfer matrix. 
We recall that the space of states of the 8VS0S model corresponds to the space of functions 
of the height s with values in the zero-weight space n[0], where H = V^^ . Note that the 
action of an operator of the type or E~^ on a zero- weight state of ?^[0] results in a new 
state of Ti which is no longer of zero-weight, and therefore does not belong to the space of 
states of the model. Hence, the only physical form factors of the model are those involving 
local operators of the type E^~^ or E^ . 

Let us first consider the matrix element of the operator E^ at site i between two Bethe 
eigenstates {{u},uju \ and |{f},a;„), where {u} = {ui,...,Un} and {v} = {vi,...,Vn} are 
solutions of Bethe equations, associated respectively to Uu and cj^. From the solution ()4.8p of 
the inverse problem, we get 

i-1 1 

= n^(^'^) • ^(^^) • n i^i^f')]''- (5-1) 

k=l k=i 

Hence we have 

// -I \jp—\f \ \ IVk'=l'^(.^k;{u},UJu) 1 ~ / N / IN 

{{u},uju\E. \{v},uj^) = —r -y > (/p^„(s)(^^„(s- 1) 

Uk=i^i^k;{v},ujy) ^ 

X (0 |CK; s-n)... C{ui;s - 1) D{^i; s) B{vi;s - 1) . . . B{vn; s-n)\0). (5.2) 



Using the commutation relation, which follows from (|2.8p . of the operator D with the operators 
B, we can express ()5.2p in terms of the partial scalar product (jS.lOp : 



{{u},uJu\E. = — f^— — > ^^Js)(/7^„(s- 1) 

\Ak=l^i^k]{v},UJ^) 



X 



Sn{{u]]{vp]i3^jU {ii};s). (5.3) 
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From (|2.14p . ()2.18p . and the value ()3.10p of the partial scalar product, we get 



Uk=i'^i^k;{u},u}u) 1 Ut=id{ut) 



{{u},uju\E- |{u},a;^) 



X 



nUi ^(^fc; {v},^v) UkKii^k - ui][vi - vk] 

n n r , 

^(-ir^+^d(.,)n \"!',i ^ E (-1)'"'^'" 

i=l '=1 ^ 5,5c{l,...,n} 

n {i-^r<v,) flint - + 1]| n fiK - - 1]| 

k^s I *=i J fce5 I " <=i J 

[L ^f^,JW\-\v\ + v,-^, + ^ + s] [s-\S\ + \S\] J 

X det [A/'-^(M;{'i)-5^^})], (5.4) 

where 7 is arbitrary, and ■0^ = if 7^ j, vj = ^j. Setting 7 = — |u| + |f | and performing 
a change of indices in the sum over s, we see that, similarly as for the computation of the 
scalar product, we can factorize out the sum over s and use the linearity of the determinant 
to recast the last sum over determinants into a single one: 

/ r„T , , I p— I r„; , , x _ Ul'Jl^iCk; {u},UJu) / 1 V- , ,.-1 h + s] \ 
{{ul,uju\hi \{v\,UJ^}-— -<j 2^ 

nfc=i^(6;W>wj l^.^L i^J J 



Uk<i i'^k - ui] [vi - vk] jri iJi - vi] 

X det rr2^({n},a;«; {■^)},a;^,)l, (5.5) 

with given by (j3.19p . Finally, the remaining sum can be seen as the expansion of the 
determinant of the sum of two matrices, one of which being of rank 1: 

{{u},u^\Er\{v},u.) = { II :rr:r:'"( ^^T >: 



T-r T{(,k;{u},uju) \ / 1 - 
nr(a,W,c..)/\L^2^^-" 

^0 



k 

Ut=i djut) 

Uk<l[Uk - Ul][vi - Vk] 



I det [Q.^{{u},Uu; {v},uj^)\ 



det [f]^({u},w„; {v},u^) - V.y{{u},uju; {v},uj^\ii)\ |, (5.6) 



where we have defined the rank 1 matrix as 



[V^{{u},Uu■,{v},u:^,\ii)\^^^ = j-' 



[7] I [Ua - ii] \Ua - + 1] 

(-i)-«(.)n{[^.-^.+i]|^^ 



(5.7) 

- <;i j^j J 

Note that we have used the Bethe equations for {v} to express V-^ as in ()5.7p . 

Remark 5.1. When {n} ^ {?;}, one can use the orthogonality of the corresponding Bethe 
eigenstates and the fact that det^ \Q.^{{u} ,u}u', {v},u}v)\ = to simplify the expression ()5.6p . 
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The computation of the matrix element {{u},uju l-E-j^^l {v},iOv) can be performed in a 
similar way, by using the solution of the quantum inverse problem for the operator in 
terms of the operator ^(^j), and by computing the action of the latter on the state | {v}, ujy ). 
One can also notice that E^"^ = 1 — E^ and use the previous result concerning the form 
factor of E^ . One obtains: 

T(^fc;{n},a;„) 1 J 1 ^-^ , .-s, . ^7 + «] 1 



Ut=i d{ut) 



Uk<l[^k - Ul][vi - Vk] 



det [0,^{{u},uju; {v},ujjj) - P^({n},a;„; {w},a;„|^i)] . (5.8) 



Combining these two results, one gets the following expression for the form factor of the 
operator af between two Bethe eigenstates: 



{{u},u}u\crl \{v},u}y) -- 
Ut=i d{ut) 



Y\k<l[^k - Ul][vi - Vk 



rdet [^^{{u},uju] {v},uj^) -2V.y{{u},uJu;{v},uJv\ii)\- (5.9) 



We recall that, in all these expressions, 7 = —\u\ + \v\. Note that these formulas are very 
similar to those obtained in the XXZ case. 



6 Master equation representation for the two-point function 

We now consider the problem of the computation of the two-point function, which can be 
tackled by summation over the corresponding form factor series. We explain here how such a 
summation can be formally performed in finite volume, hence resulting into a multiple integral 
representation for the two-point function, similar to the representation that was introduced 
in |60j in the XXZ case under the name of master equation representation. This terminology 
was due to the fact that such a formula could be understood as the common result of two 
possible ways of computing the two-point function as a sum over elementary objects: over 
form factors on the one hand, or over elementary building blocks on the other hand. 

We consider here a quantity which is the analog of the generating function studied in [60] . 
Namely, for a generic complex number k, we introduce the operator 



= n (^ + ^"i) = n {^=r + -•^ry 



(6.1) 



The mean value of this operator in the ground state | {u},u}u ) of the transfer matrix, 



{{u},UJu\{u},UJu) 

which can be seen as a polynomial of degree m in k, is a generating function for the two-point 
function of the finite-size one-dimensional quantum dynamical model whose Hamiltonian is 
obtained from the finite-size transfer matrix. In the thermodynamic limit, it also gives the 
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probability, in the CSOS model, that two sites on a same Hne at distance m have a difference 
of height £ < m: the latter is given by the coefficient of of ()6.2fi 

The quantity (|6.2p can be evaluated by summation over form factors corresponding to a 
suitable basis of the space of states. In the present case, such states can be conveniently chosen 
as the eigenstates of the k- twisted transfer matrix t^iu) = A{u) + k,D{u). The latter can be 
constructed similarly as in the untwisted states of the form (l2T3]l - (fmD or (fZTHD - 

(|2.19p with parameters {v}^ satisfying the following system of K-twisted Bethe equations: 

[.^_,. =(-1) <^^-m^—^^ J = l,...n. (6.3) 

The corresponding eigenvalues are 

Tk{u;{v\,uj)=uj a{u)\\— ^+ -1) Ku d{u)\\— — . 6.4) 



[VI - u\ fj^ [u - Vl\ 



Indeed, from the solution ()4.9p of the inverse problem, the operator ()6.ip can be expressed 
as follows: 

m 1 
i=l i=m 

Hence, provided the eigenstates \{v}^,ujv) of t,i{u) form a complete basis of the space of 
states, one can expand the mean value (|6.2p as 

lO^ \ = V (i^i'^" K^}-^'^^) \Ql,m\{V'}^^u) 

^ -TT Tk{^u {v}^,'^v) {{u},^u I {v}^,i^v) {{v}^,^v I {^},^n) ,g 

ri -'■i T"te;{^*},Wn) {{u},^u\{u},U}u) {{v}^,U},\{v}^,UJ,)' 

and then use the determinant representations obtained in Section [3] to express the scalar 
products. 

In the untwisted case, for a rational parameter 1] — z ^^'^ ^'^^ ^ ~ -^/^' completeness 
of quasi L-periodic Bethe eigenstates ()2.13p - ()2.14p with generic multiplier a (i.e. associated to 
io such that (— l)"w^ = a) and generic inhomogeneity parameters {^} was shown by Felder, 
Tarasov and Varchenko in [27], at least for any large enough odd integer L > n. Their proof 
can easily be extended to the case of L-periodic K-twisted Bethe eigenstates with arbitrary 
K (at least in a vicinity of k = 0, which is enough for the study of the polynomial quantity 
(j6.2p ). and to more general rational parameters rj = r/L for L odd and large enough. This 
problem is discussed in Appendix [El Hence, in (j6.6p . the sum should be understood over all 
admissibl^ off-diagonaH solutions {v}^ of the K-twisted Bethe equations ()6.3p with n = N/2 
(N = 0), associated to all possible values of uj, £ C such that {—1Y"'uj^ = 1. 



®We suppose here that m <^ L. In fact, in this section, we are not interested in the specific properties 
of the root of unity case, but rather in obtaining a formula that could be extended by continuity to the case 
of irrational 77. Therefore, when necessary, we may suppose L large enough (for instance L > N) to avoid 
problems of cyclicity on the lattice. 

solution {v} of the system H6.3|) is called admissible, if Hfe^i l^i ~ Cfc] YTi=i — f ! + 1] 7^ 0, j = 1, n. 

'^See footnote 111 
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The computation of the scalar product performed in Section [3] can easily be generalized 
to the K-twisted case. Using the same arguments as for (|3.18|) . one obtains that the scalar 
product between a k- twisted Bethe eigenstate ( {v}^ , I and a general state | {w}, tOw ) of the 
form (|2.13|) - (|2.14p with n = N/2 can be represented as|j 



< b + s] 



I H J Wj<k['>^j-Vk][Wk-Wj\ 



X dei[^i^\{v\,io,-{w},u^% (6.7) 



with 7 = — |f I + \'w\, and 



1 ][Vi-Wj+-i] UJ^[Vi-Wj+-i+l\\ -A- 



[7] bi-^^j] 



[Vi - Wj + 



_^ J_ j - Wj + 7] cj^ [uj - Wj +j - 1][ K 



[7] I [^i - Wj] Uw [Vi - Wj 



[vt-Wj - I]. (6.5 



The "square of the norm" of a K-twisted Bethe eigenstate is given by (j3.20p or, alternatively, 
in terms of the function 



3^k;w('u; {v}) = a{u) '^[vi - u + I] + Ku ^ d{u) '^[vi - u - I], 



(6.9) 



as 



1 ntidjvt) 

{-P\ri[j^k[v3-Vk] 



det 



dvL 



yK-L,.Avj\{v}^] 



(6.10) 



Hence, using these determinant representations for the scalar products, one can rewrite 
(EH) as 



[7i> + s] 



T 



E 



detn [^ly^{{u},uJu;{v}^,uJv)] • det„ [rir^J{7;}^,a;„; {n},a;„)] 



, (6.11) 



in which we have set 7^, = |f | — |n|. As in [60j, the sum over admissible off-diagonal solutions 
of the K-twisted Bethe equations in (|6.11|) can be reexpressed as a multiple integral, and one 
obtains the following representation (master equation) for the generating function ()6.2p : 



n 



d'^z -Q r«(^i;{z},a;) 



r(MJ 



Ew^ [iz + s] 



det„ [Q^^({n},a;„; {2;},a;)] • det„ [Vr\^({z},LO]{u},uJu)\ 



(6.12) 



Here we only give the case n = N/2 (i.e. K — 0). In the general case, there is an additional sign in the 
expression of o!"^^ as in (|3.19|l . 
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The sum in ()6.12p is taken over all u; € C such that {—lY"'uj^ = 1, and the integration contour 
is such that it surrounds (with index 1) all poles corresponding to solutions {v}^ with n = N/2 
of the K- twisted Bethe equations associated to uj. The other poles of the integrand, which lie 
outside the contour, are poles at Z£ = i = 1, . . . ,m (issued from the factors TK{^i; {z},uj)), 
and poles at zi = Uj, j = 1, . . . ,n (contained in det^ ['^^-^^ f^u)] )• By considering 

these poles, one can, as in the XXZ case |6Dj, rewrite the above integral in the alternative 
form 

/o- \ = lM!l!v / j^frT^^(e^^{£M / >p u;^ bz + s] Y 



"r({c})ur(W) '-^ecL 



det„ [Q^^{{u},Uu;{z},u;)] • det„ [Q^/^}^^{{z},uj; {u},uju 

detn [Qo{{u},Uu;{u},UJu)] ■lYj = iyK;uiiZj;{z}) 



(6.13) 



7 Conclusion 



In this paper we have generalized the algebraic Bethe Ansatz approach to correlation functions 
to the dynamical case, i.e. to the case of an i?-matrix depending on a dynamical parameter and 
satisfying the dynamical Yang-Baxter equation. In this framework, we have obtained finite- 
size representations for scalar products of Bethe states, form factors of local operators and 
two-point correlation function. These results hence initiate the study of correlation functions 
of the (cyclic) 8VS0S model within the framework of algebraic Bethe Ansatz. Applications 
of these formulas at the thermodynamic limit will be considered in a further publication. 

As we have seen, one of the difficulty of the ABA approach to correlation functions is 
that it relies on the existence of a compact and manageable expression, preferably as a single 
determinant, for the scalar products of Bethe states: in practice, such a representation may 
not be easy to obtain, and it is even possible that it does not exist for all models, at least 
in the simplest expected form of a single determinant (see for instance the recent preprints 
[21 [8] concerning S'f/(3)-invariant models). In the case of general SOS models, in fact, it seems 
that we can only obtain representations of height-dependent (partial) scalar products (j3.2|) as 
sums of determinants. As in the case of the partition function with domain wall boundary 
conditions |84[ I81j. this is due to the shifts undergone by the dynamical parameter. However, 
as mentioned in Section [3l it is not very surprising that shifts of the dynamical parameter 
prevent us from representing the height-dependent partial scalar product (and hence the 
partition function) as a single determinant, since such shifts also prevent the corresponding 
states, i.e. states of T-L[0] of the form (j3.15p . from directly being Bethe states. In fact, when 
one considers true Bethe states as functions of the dynamical parameter, and their natural 
scalar products which, in the cyclic case, can be expressed as a finite sum over all values of 
the dynamical parameter, one obtains for the latter a single determinant representation very 
similar to its XXZ analog. It follows in particular that the norm of a Bethe eigenstate can be 
expressed as a single determinant. As we have seen, this opens the way to the computation 
of form factors and correlation functions. 

As a final remark, we would like to mention an interesting and promising alternative 
method: Sklyanin's separation of variables (SOV) method [881 189j . which has been recently 
developed in the direction of the computation of form factors |41[ [80l I79j . The advantage 
of this approach is that it enables one to free oneself from the main difficulty of the ABA 
approach to correlation functions, namely the fact that, depending on the model we consider, 
we have a priori no guaranty that a convenient expression for the scalar products does exist 



20 



and — should it exist — no model- independent constructive method to derive it. In fact, 
in the SOV approach, it seems that scalar products of separates states can automatically be 
expressed as a determinant. Let us note however that, at least at the present stage of the 
art, the determinant representation obtained from ABA still presents some advantages, the 
main one being probably that it is more appropriate than its SOV counterpart for taking the 
thermodynamic limit. 
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A Definition and useful properties of the function u [u] 

Let us set p = e^'^*'^, with > 0. Let 6i{z;t) be the Theta function of quasi-periods 1 and 
r and nome p^^"^ = e*'^'^: 

oo 

ei{z; t) = 2p^/^ sin TTz (1 - 2p" cos 2ttz + _ p") 

n=l 

oo 



-I 

n=—oo 



In the generic SOS model (i.e. for -q generic), the elements (|2.3p - (j2.4p of the dynamical 
i2-matrix are usually defined in terms of the function 

[n] = 0i(7?n;T), (A.2) 
which is odd, entire, and satisfies the quasi-periodicity properties 

[n + 1/r/] = - [u] , [u + t/t]] = -g-^^^e-^*^''" [u] . (A.3) 
This function also satisfies the following useful properties: 

1. Addition formula: 

[x + u][x — u] [y + v][y — v] — [x + v] [x — v][y + u] [y — u] = [x + y] [x — y][u + v] [u — v]. (A. 4) 

2. Let t G C and / be of the form 

/(x) = ygo-) ^'^""V|'''^^"°^+-^ , (A.5) 

r [x-b?]...[x-b?] 

where ap'\ 6^, C7(j) G C such that a[^^ + ■■■ + a^/}^ - b^f b^P = t mod 1/??, Vj. 

Then / satisfies the quasi-periodicity properties 

/(X + 1/r/) = (-1)'"/(X), f{x + r/r?) = ^-l)m ^2.^r,t-i.m{2r,.+r) ^ (^.6) 



21 



Moreover, if / is entire (i.e. if the singularities at x = b^p are all removable), then it 
can be factored as 

f{x) = C[x-ai]...[x-am], (A.7) 
where C, ai , . . . , am are complex numbers such that ai + ■ ■ ■ + am = t. 

3. Frobenius' determinant evaluation: 

det (^zyi±^\ = - \y\ + 1] "^^-</i"'";" " , (a.s) 

l<i,3<n\ [Xi-Vj] J Ulj=i[xi-yj] 
where |x| = xi H h x„, \y\ = yi -\ h yn- 

4. For any positive integer L, one has the following identity: 

where [u]^ = 9i{rju; Lt). (|A.9p can be shown by relabeling the sum {k k + Lj) in the 
Fourier series expansion: 

fe=— 00 

In the cyclic case, when Lr] = ri + with L,ri,r2 integers, it is convenient, following 
Baxter [6i^j, to re-define the function u 1-^ [u] as 

[n] = 0i(??u;T)e^""2i"'. (A. 11) 

The modified function (jA.lip is odd, entire, and satisfies the quasi-periodicity properties: 

[u + 1/rj] = -e'"'?(2"+^) [u], [u + T/r,] = -e~*"*(2«+|) j^^^ ^^_^2) 
It is also quasi-periodic of period L with a mere sign: 

[u + L] = (-l)'-i+^2+'^i'-2 [n], (A.13) 

so that a ratio of two such functions is periodic of period L. The modified function (|A.lip 
still satisfies the properties 1. and 3. above. If moreover / is of the form (jA.Sp with a^"'^ + 
■ ■ ■ + fll+m ~ ^1"''* — ■ ■ ■ — b'^P = t mod L, then / satisfies the quasi-periodicity properties 

f(x + 1/r?) = (_l)-e*-*(2— 2*+-^)/(a;), /(a; + r/r^) = (-l)-e'"*(''""-'*+"^)/(x), 

and can be factored as in (jA.7p . In that case, it is enough to find m independent zeroes to 
prove that / vanishes identically. 



^"Equation (9) of paper I of [6] corresponds in fact to L77 = 2ri +2r2r. Here we get rid of factor 2 by allowing 
a sign in the L-periodicity (|A.13|l of the function it 1— !> [it] . 
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In fact, as noticed in [23], the modified function (|A.lip is nothing else but a usual theta 
function of the form (jA.2p with different quasi-periods. This can easily be seen from the 
transformation property of the theta function (|A.ip under the action of the modular group: 

e-'^.e,(..r) oc f„. aa.y (° ^) e SL(2,Z). (A.14) 

Hence the modified function (jA.lip is simply proportional to 

[u] cx 6i{r]'u;T'), with rj = — and t' = ^— . (A.15) 

L ri + r2T 

In this expression, r is the greatest common divisor of r\ and r2, and the integers a,b,ri,r2 
are such that ri = rri, r2 = rr2 and ari — br2 = 1. Hence all the properties of ()A.2p can 
trivially be extended to (jA.lip with the mere replacement of rj by r]' and r by r', and the 
study of the case Lij = ri + r2T is therefore equivalent to the simplest rj = r/L case. 



B Partition function of the 8VSOS model with domain wall 
boundary conditions 

In this appendix, we recall the result of [83] concerning the representation of the partition 
function of the 8VS0S model with domain wall boundary conditions as a sum of determinants, 
useful for the computation of scalar products. For a inhomogeneous square lattice of size NxN 
with parameters ui, . . . ,U]\f,^i, . . . ,^Ar, this partition function can be expressed in terms of 
the entries of the monodromy matrix ()2.7p (with inhomogeneity parameters ^i, . . . ,£,n) as 

ZNiui, ...,un;Cu...,(n;s) = {0\C{un; s) C{un-i;s + 1) . . . s + (A^- 1))| ), (B.l) 

where | ) = e_ (8> • • • (8) e_. It can easily be shown that this quantity is uniquely determined 
by the following properties \84:\ [HT] : 

2. Zj\f{ui, . . . ,uj\f', Ci; • • • ) Cn', s) is symmetric in ui, . . . ,uj\f and in ^i, . . . , ^tv separately; 

3- n^^=i[^a ~ + 1] Zj^{ui, . . . ,uj\[;^i, . . . s) is a holomorphic function of each Ua 
and each 

4. it satisfies the recursion relation: 

Zjv(ni, . . . , txiv; 6, ■ ■ ■ ,^jv; s)L^=^^_i = + ^ _ 1] 

^ T t r~rr -^Zjv_i(ui,...,nAr_i;Ci,...,CiV-i;s); (B.2) 

- Cjv + 1] \ik - iN\ 

note that we have also 

ZAr('Ui, . . . ,'UAr;^i,. . • , Ca^ ; = ZAr_i(ni,. . . ,UAr_i;^i, . . . ,CAr_i;s + l); (B.3) 
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5. Zn = /(til) as a function of ui satisfies the quasi-periodicity properties: 

/(m + 1/??) = e-2-¥(^+^)/(ni), fim + r/r/) = e-2'''^x(-+^)/(ni). (B.4) 

In |84^ [ST] , a representation of the partition function as a sum over the permutations of the 
set {^} was proposed. In [8l], it was shown that this representation could be reduced to the 
following sum of determinants: 

Zn[{u};{^};s] - 



[in| - i^j + 7 + s + iv] n,<fcK- - ^fc] [Cfc - C,] 

X E (-1)'^' ^^tltr-rJlf'^ detK({n + c^^};m)]. (B.5) 

In this expression, 7 is a generic parameter (the left-hand side being independent of 7), Af^ 
is given by ()3.1ip . \u\ = ui + ■ ■ ■ + un, |CI = ?i + ' ' ' + ^N, the sum runs over all subsets S of 
{1, . . . , A^}, with IS"! being the cardinality of 5 and {u + 5^} = {uj — 5|'}i<j<Ar with 

6^ = 1' '''l'^ (B.6) 
' \0 if J ^5. 

Formula (|B.5p was proved in [M] by induction over A^. It is easy to see that, following the 
lines of the proof of [H], one can also obtain the equivalent formula: 

ZN{{u};{i}]s) - 



X E (-1)'^' ^^tltr-Jf'^ detK(M;{C-^n)]- (B.7) 

Remark B.l. The partition function ZAr({M}; {^}; s) has no pole at \u\ — \^\ + 7 + s + = (it 
is holomorphic except at the points Uq, — ^^ + 1 = 0), which means that the sum of determinants 
in the second line of (jB.Sp or (|B.7p has a zero at this point. 
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C Solution of the recursion relation ( 13.81 ) for the scalar product 

We show here that the solution of the recursion relation (|3.8|) . (|3.9|) . with the initial condition 
(|3.7p . is given by the following expression: 



Gi%,,.../„{{u}; {vi, . . . ,vk}] s) 



[s-k] 



k-1 

n [s + iv - 2n + j] 

j=o I- ■'^ l=k+l 



n 



[7]" [H - \S,i\n-k - \v\k +1 + s] 

1 n-=iK--c£. + i] 



TT i^ii - + 1] A 



11^=1 4^j)n"=iQ'(^i) Wk<j<i[^tj 



3^1 

( 1^51 + 151 [7 + ^-|g| + |5|] 

sctJ [^-\S\ + \S\] 

Sc{l,...,k} 



ie5 



1=1 



det [J\f^{{u}; Vi 



ss 



SS t-S 



5 ^fc ) ?^ 



fc+1 ' 



(C.l) 



S5 



^? -^f^' =?^, -'^f with 



Here we have used the notations of Proposition 13. II with v 

Sf given by (|B.6p. = vi -\ \-Vk, \£,e\n-k = 6+1 H ^ Cn- 

Let us first remark that (jC.ip provides the correct expression for G^^^ ^ : this follows from 
(|3.7|) and from the representation (|B.7p for the partition function with domain wall boundary 
conditions. 

Supposing that (jC.ip is valid at level /c — 1, we now want to prove that it is also valid at 
level k. From the recursion relation (13.81). (13.91). we obtain 



Gi%,,...,e„{{u}; {vi, . . . ,vk}; s) 



n 



[s - j] 



1 



n 

l=k+l 



[^i, - + 1] i^i, -vk-i] u]=i K - + 1] 



N 



lYj=k+i[Ce, - Ce, + 1] 



n ■=! rf(^i) n"=i ^(^.) n.<i<j6. - 

Ui<iWj - ui]Uj<i<k[vi - vj] 



N 

n 



\ [vk - Cj] 



■f^'Hvk), (C.2) 
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with 



k 



ni=iK-e4 + l] / .^|5|+|5| [7 + g-|'g[ + l^|] 



E (-1)' 



5c{l,...,fc-l} 

X n ni". - + 111 n ((-D'v-ni- - "> - n) 

j=i L /=i J \ ) 

X det [M^{{u}; vf^..., vl^, , ef,, • • • , Cl )] • (C.3) 

Note that the sum in (|C.3p (which corresponds to the sum in ()3.8p ) has been extended to 
all values of ik- indeed the sum over determinants, which is antisymmetric by exchange of 
two different vanishes for = ^j; j = k + 1, . . . ,n. The function (|C.3p of is an entire 
function of the form ()A.5p with m = N and t = ^^Li — s — N + 2n — k. It can also be 
written as 

f^'Hv,) - TTf.. - £-1 V [1] [. + iV-2n + fc + .,-6J fr [6. - + 1] 

1 E (-1)'"^'" 



ni=fc+lK^, - ?4 - 1] • U^li^ik - Vj\ 5c{l,...,n}\{fe} 

Sell,.-,*:-!} 

X n s + 1] n (-d^-^'^ m-^ - - n 

j=i y 1=1 ) .^g I 1=1 



X det [M\^^ {{u};vf^..., vl^, , £4 , ^g^^ , . . . , Cfj] , (C.4) 
with [ACj''^] = [My]^j for j / A: and 

n 

i=i 

[7 + 5-|5| + |g|] [n,-£4.+7] _ [7 + .-|5| + |g|-l] - £4 + 7 + 1] ) 
+ [ui-Ce,] [s - \S\ + \S\ - 1] [n, - £4 + 1] / 

TTr [1][7][^. - £4 + s - |5| + |S| + 7]K - £4 - . + |5[ - |g| + 1] 

f=\ ' " [s-\S\ + \~S\\[s-\S\ + \~S\-l][u.-i,,] 
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in which we have used ()A.4p . This function is equal to 



S'C{l,...,n}\{fe} 3=1 I ^ ■''^ /=1 ) jeS ^ '=1 

5c{l,...,fc-l} j^S 



N 



X 



- + 1] det [AC;^)({n}; ^ . . . , vl^„Vk, ^g^^ ^fj] 



iV 



i=i 



l.rK -2TTr _ . l V- [s + N - 2n + Vk - ^j^ + k + 1] 1 



i=i <j=fc+i 

n ^ k—l ^ k—l ( / \ ^ 

j=fc+l L'^^^ S,«5J j.^^ LS«, J ^ 5C{l,...,n}\{fc} i = l I ^ •'M=l 



•J'^? sc{i,...,fc-i} 



X 



n (-^y'^u' U^ni - V, - 1] det [AAf ({«}; ^ . . . , ^f^^, - 1, C,, • • • 



y. [S + N -2n + Vk-Vq + k] Uf=l[Vk - ^j] yj 1 -A- 1 

^ [hi - IC^In-fc - l^-U-l -'(^9 + 7 + 'S] [Wfc-I'g] f-Jl [Vq - Vl] [^e^ - Vq - 1] 

E n HN - V, + 4 n I (-1)*-.-= Ill'- - v> - 111 



J6S 



/ = 1 



5c{l,...,n}\{/c}i=l I ^ ^' 1=1 
5c{l,...,fc-l} jY-S 

- (-l)'-Vdet [AA«(M;^f ^ . . . ,vl!,,Vq,il^^, . . . |, (C.5) 

with [T^f (M;.f^...,.f^,,.,,ef,,^,...,efj]^^. = K(M;.f^...,.fA,et...,efj],, 
for j ^ k and 

[ATf (M;.f^...,.fA,..,et,,...,efj],, 

= _-rr[^ I] [^][7][n» - T^fc + 5 - |5| + |5| + - - , + |5| - |S| - 1] 
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fru _^ I] f h + ^-\s\ + \s\] [u^-vk + l] 

fJl ' ' ^ I - |5| + \~S\\ [ui - vk] 



[^ + s-\S\ + \S\ + l] [u,-vk + i-l] 
[s-\S\ + \S\ + l] [u^-vk-l] 



, (C.6) 



in which we have used ()A.4p . Note that, due to the fact that {u} satisfies the Bethe equations, 
there is no pole at Vk = Ui in (jC.Sp . Note also that, due to Remark (jB.ip . there is no pole 
at \u\ — \S,i\n-k ~ \v\k + 7 + s = since the corresponding sums over determinants have a 
zero at this point. Therefore the function (jC.Sp is an entire function of the form (jA.Sp with 
m = N and t = XljLi Cj ~ ^ — N + 2n — k, and it takes the same value than (jC.4p at the 
independent points Vk = ^j, j = I, ■ ■ ■ , N. Hence the equality between (|C.4p and (|C.5p . 
Up to the product nj=i[^fc ~ ^il' second term in ()C.5p is equal to 

(,yH ^2 V [s + N-2n + Vk-^i^+k + l] nr=lh-g^J TT 1 



k-l , 



X y (-i)i^ili±i^^^iJfclldetrAL(M;t-f^...,^f^i,ef,e£ ,...,ef )], 

which is equal to zero since the sum over S is antisymmetric with respect to the variables 
and since the variable appears 
the third term in (jC.5p is equal to 

fc-l r , i^r ^ . ,71 - k-l 



and since the variable appears twice in this sum. Also, up to the product nj^ibfe ~ Cj] 



^ [s + N - 2n + Vk - Vq + k] 1 -j-j- 1 y-j- 1 

^ IH - \^e\n-k - \v\k-l -Vq + -f + s][Vk- Vq] [Vq - Vl] ^. j-^^ ^ [^e, - Vq - I] 



X E n |^n[«'-"i+4n|(-irwn['''-".-4 

5c{l,...,n}U{g}\{fc} je{l,...,fe-l}l ^ ■'^=1 ) j&S^ l=i ) 

Sc{l,...,fc-l}U{g} U{g}\S 

which is equal to zero due to the antisymmetry of the sum over determinants by exchange of 
two (the variable Vq appears twice in this sum). In the two above expressions, the notation 
{1, . . . , n} U {q} or {1, . . . ,k — 1} U {q} means a set that contains twice the element q. 
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Finally, we obtain, 

,(fc). . ^ [s + N-2n + k] Uf=i[vk-(j] 

^"""^ " M - ledn-fc - \v\k +1 + S] U]=k+l[^i, -Vk-l]- Uj=l[vk - V,] 

5c{l,...,n} J = l I ^ ■'M=l J j^s ^ «=1 ^ 

5c{l,...,fe} j^S 

which means that Gf^J^^ ^^({n}; {vi, . . . , Vfc}; s) is given by (jC.ip . 

D Orthogonality of two different Bethe eigenstates 

We show in this Appendix that, similarly as what happens for the XXZ chain [59j, the 
matrix Q^{{u},uju', {v},ujv) (|3.19p giving the scalar product of two different Bethe eigenstates 
{{u},uju I and | {v},ujv ), with {u} and {v} two different off-diagonal solutions of the Bethe 
equations (|2.15p . admits a zero eigenvalue. 

Let us first suppose that Uj ^ Vk,yj,k. Then, using the Bethe equations for {v} and 
factorizing the product {—^Y^a{vj)Y[]^=i[ui — Vj + 1] out of each column j of the matrix 
Q.y{{u},uju', {v},ujy), we obtain a matrix M^{{u},uju', {v},ujv) of elements 

[7] \ [ui - Vj] [ui - Vj - 1] j ujI [vi - Vj - 1] [ui -Vj + iy 

The action of the transpose of the above matrix on the non-zero vector w of components 

[h=i[ui-ui\ 

produces a new vector w of components 

n 

Wj = E [M^{{u},iOu;{v},UJvY]^.^j 

J I Wi) J fj^ [vi -Vi-l\ [Ul -Vi + l\ 



with 



^ [7] [tx, - + e] nr=iK-«d 

fJl [vi -ui-e\ 
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The equality between (|D.4p and (|D.5p follows from the residue theorem applied to the elliptic 
function 



n 



[z-Vi+-f + e\ [z - vi\ 

integrated on an elementary cell (we recall here that 7 = \v\ — \u\, which means that the 
function (|D.6p is effectively doubly periodic, and that its integral on an elementary cell 
vanishes). Replacing the value (jP.Sp of into the expression ()D.3p . we get that w = 0, 
which means that the non-zero vector w is an eigenvector of M* with zero-eigenvalue, and 
hence that the scalar product ( {u},uJu \ {v},uJv ) vanishes. 

In the case when some (not all) of the variables uj coincide with some of the variables Vk, 
the proof is similar but one has first to take the limit in the corresponding columns of the 
determinant. 



E Completeness of Bethe eigenstates 

In the untwisted case and for a rational parameter V — the completeness of Bethe eigen- 
states (with arbitrary multiplier through the transformation s — >■ s + L) was studied by 
Felder, Tarasov and Varchenko in [27j . More precisely, there was shown the following result 
(Corollary 19, Theorem 21 of [27]): 



Theorem E.l. 127}/ Let N = 2n and rj = 1/L for some odd integer L > n. For generic a, 
{^}, there are d = Ldim?^[0] admissihl^^ off-diagona^ solutions of the Bethe equations 

with uj^ = (— l)"a such that the corresponding Bethe eigenstates ()2.13p - ()2.14p form a basis of 
the d-dimensional vector space of quasi L-periodic functions f G Fun(?^[0]) with multiplier a, 
i.e. such that f{s + L) = af{s). 

For our purpose, we need a slightly modified version of Theorem IE. II 

Theorem E.2. Let N = 2n and ij = r/L for some relatively prime integers L,r with L 
odd and L > n. There exist kq > such that, for < \k,\ < kq and generic {^}, there are 
d = Ldim?^[0] admissible off-diagonal solutions of the n-twisted Bethe equations 

n n 1^^^^^ = .c.-, j = l,...,n, (E.2) 

with oj^ = (—1)''"', and the corresponding K-twisted Bethe eigenstates (|2.13p - (j2.14p form a 
basis of the d-dimensional vector space of functions f G Fun('?^[0]) such that f{s + L) = f{s). 

"See footnote El 
^^See footnoted 
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The proof of Theorem lE.ll in [27] follows from the study the following eigenvalue problem 
with multiplier condition (this problem is related to the study of solutions of qKZB equations 
for the elliptic quantum group Er^n{sh))- 

i7,(e)V' = e,>, j = l,...,N, 
^peFun{n[0]) such that ^p{s + L) = aip^s). 

Here Hj(^) are commuting difference operators defined in terms of i?- matrices of the elliptic 
quantum group Er^r]{sl2)'- 

^jiO = RjJ-i ■ ■ ■ Rj,iTs^ Rj,N ■ ■ ■ (E-3) 

where ^ = (^i, . . . , ^^r) is a fixed generic point in C-^, and where we have used the shorthand 
notation 

fc-i 

Rj,k = Rj,k {^j -^k;s + Y,hi). (E.4) 

1=1 

In particular, in the fundamental case, the difference operators (|E.3p coincide with special 
values of the transfer matrix of the SOS model, and their eigenfunctions (for generic 
coincide with Bethe eigenstates. 

The first part of the proof of [27] consists in constructing common eigenfunctions to the 
operators Hj{S^), j = 1, . . . ,n. To this aim, the spaces dual to tensor products of evaluation 
Verma modules over ET-^r){sl2) are realized as particular spaces of meromorphic functions of 
n + 1 complex variables (typically the n spectral parameters and the dynamical parameter) 
with adequate quasi-periodicity properties. A basis of this space of functions is explicitly 
identified. Common eigenfunctions (with given multiplier a) of the commuting difference 
operators Hj{^) are then constructed in terms of these basis elements evaluated at the solutions 
of the system of Bethe equations (|E.ip . These eigenfunctions are identified (up to a factor 
that vanishes for diagonal solutions of (jE.ip ) with Bethe eigenstates. We refer the reader to 
[27] for explicit details about this construction, which is valid for any value of rj and of the 
multiplier a. 

The second part of the proof of [27J consists in showing that, in the case n = A^/2, 
T] = 1/L, and for generic a, the previously constructed set of eigenfunctions is indeed complete. 
For a large but finite, one can identify Ldim?^[0] different sets of admissible off-diagonal 
solutions of the system (|E.ip . Using the explicit expressions, previously obtained through 
the aforementioned construction, of the corresponding eigenfunctions, the authors of |27j are 
finally able to show that the latter are indeed linearly independent, at least for a large enough 
(see [27] for details). 

To prove Theorem IE. 21 one needs to extend these arguments to the K-twisted case with 
multiplier 1. It is in fact easy to see that, in the previous reasoning, one can use k as a 
free parameter instead of a. We use the fact that solutions of the system of Bethe equations 
at K = are very simple and that, for L odd, there exist L possible distinct values of 
the parameter uP' such that oj^ = (— l)*^". Hence, considering the system ()E.2p when k 
tends to zero and using the implicit function Theorem, one can identify Ldim?^[0] differentia 
sets of admissible off-diagonal solutions {v^j'^{n)}i<j<n (corresponding to all possible subsets 

'^^up to periodicity in I/77 and t/t). The condition L odd and L > n ensures that the zero weight space "^[0] 
reduces to the space corresponding to n = N/2. 
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I = {ii, . . . C {1, . . . , N} and all possible integers i € {0, . . . L — 1}) which are continuous 
functions of k in a vicinity of k = 0. They are such that 

v^/{k) = - 1 + Ke^'^r n}(0) + o(k), (E.5) 

with 

To show that the corresponding Bethe eigenstate form a basis of the corresponding space 
of functions, i.e. that they are linearly independent, one can adapt the arguments of [27J to the 
K-twisted case by explicitly constructing common eigenfunctions to the K-twisted commuting 
difference operators 

H^;\0 = Rj,-i---Rj,ir^'K^,^^Rj,N...Rj,j+i, with K(-^=(l^ J!). (E.7) 

The linear independence, for k. small enough, of the eigenfunctions corresponding to solutions 
(jE.Sp then follows from the same reasoning as in [27j, which can easily be extended to the 
case of rational rj = r/L. 

Alternatively, one can remark that: 

• the scalar product (|6.7p of two Bethe states corresponding to different off-diagonal solu- 
tions of the K-twisted Bethe equations ()E.2p vanishes (the proof of Appendix|D]is indeed 
directly generalizable to the k- twisted case). 

• the "square of the norm" (|6.10p of a K-twisted Bethe state corresponding to an admissible 
off-diagonal solution of ()E.2p is non-zero in a vicinity of k = 0. Namely, one has, 

for {v] = {w^'^} given by (iKSj) with co = e^'^'s'-^^'^r and Uj{0) = lij(O) given by ^Elh . 

This prove the linear independence, for < |k[ < kq for some kq > 0, of the Ldim'H[0] 
K-twisted Bethe eigenvectors corresponding to the solutions (jE.Sp . and hence ends the proof 
of Theorem IE.2I 
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